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Hub-dominant  matrices  are  natural  extensions  of  hub  matrices.  In  this 
article  we  study  eigengaps  of  the  Gram  matrix  associated  with  a 
hub-dominant  matrix.  A  class  of  hub-dominant  matrices  is  then 
constructed  by  using  equiangular  tight  frames. 
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1.  Introduction 

The  notion  of  hub  matrices  was  first  proposed  by  Kung  and  Suter  in  [1].  A  matrix  is 
called  a  hub  matrix  if,  without  loss  of  generality,  the  Euclidean  norm  of  its  last 
column  (called  hub  column)  is  greater  than  that  of  the  rest  columns  (non-hub 
columns)  and,  in  addition,  all  non-hub  columns  are  orthogonal  to  each  other  with 
respect  to  the  Euclidean  inner  product.  The  framework  of  hub  matrices  can  be  used 
to  describe  a  variety  of  wireless  communication  systems.  For  example,  the  hub 
matrix  theory  was  applied  to  beamforming  MIMO  communication  systems  in  [1], 
The  eigenstructure  of  the  Gram  matrices  of  the  corresponding  hub  matrices  was 
exploited  to  analyse  the  performance  and  capacity  of  the  systems.  Here,  the 
eigenstructure  of  a  positive  semi-definite  symmetric  matrix  refers  to  the  ratio  of  its 
two  leading  eigenvalues.  One  mathematical  issue  in  hub  matrix  theory  is  to  infer  the 
eigenstructure  of  the  Gram  matrix  of  a  given  hub  matrix  through  the  Euclidean 
norms  of  all  columns  of  the  hub  matrix  and  the  Euclidean  inner  products  between 
the  non-hub  columns  and  the  hub  column.  For  instance,  the  eigenstructure  of  the 
Gram  matrix  of  a  hub  matrix  was  characterized  by  the  hub  column  and  a  non-hub 
column  with  the  largest  Euclidean  norm  in  [1].  In  our  recent  work  [2],  an  improved 
estimate  on  the  eigenstructure  has  been  achieved  by  fully  exploiting  information 
from  the  hub  column  and  all  non-hub  columns. 

The  Gram  matrix  of  a  hub  matrix  is  a  so-called  arrowhead  matrix  whose  all 
principal  minors  except  itself  are  diagonal  matrices.  In  other  words,  the  Gram 
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matrix  can  be  obtained  by  bordering  a  diagonal  matrix  with  one  column  from 
the  bottom  and  one  row  from  the  right.  This  diagonal  matrix  can  be  extended 
to  other  types  of  matrices  such  as  block  diagonal  matrices  and  diagonal 
dominant  matrices.  Actually,  when  the  diagonal  matrix  is  replaced  by  a  diagonal 
dominant  matrix,  the  corresponding  hub  matrix  becomes  a  hub-dominant  matrix 
(see  [1]).  Hub-dominant  matrices  can  be  used  in  areas  such  as  distributed 
beamforming  and  power  control  in  wireless  ad  hoc  networks  [1].  The 
eigenstructure  of  the  Gram  matrix  of  a  hub-dominant  matrix  has  been  studied 
in  [1].  In  this  article,  we  will  improve  the  estimate  of  the  eigenstructure  based  on 
our  recent  work  [2].  We  will  also  construct  and  analyse  a  class  of  hub-dominant 
matrices  via  equiangular  tight  frames. 

This  article  is  organized  as  follows.  In  Section  2,  we  give  the  definition  of 
hub-dominant  matrices  and  review  the  result  on  the  lower  and  upper  bounds  of 
the  Gram  matrix  of  a  hub-dominant  matrix  in  [1].  Few  remarks  on  the  result 
are  made  to  motivate  our  current  work.  In  Section  3  we  develop  some  useful 
lemmas  for  estimating  eigenvalues  of  positive  semi-definite  symmetric  matrices. 
In  Section  4,  we  present  our  new  results  on  the  lower  and  upper  bounds  for  the 
eigengap  of  the  Gram  matrix  of  a  hub-dominant  matrix.  The  improvement  of 
the  bounds  over  the  result  given  in  [1]  is  discussed.  In  Section  5,  we  construct 
a  class  of  hub-dominant  matrices  through  equiangular  tight  frames.  Our 
conclusion  is  drawn  in  Section  6. 


2.  Hub  and  hub-dominant  matrices 

In  the  original  work  of  [1],  the  definition  of  hub  matrices  was  introduced  first,  and 
then  extended  to  a  general  case  of  hub-dominant  matrices.  The  present  presentation 
differs  in  that.  We  first  give  the  definition  of  hub-dominant  matrices  and  then  treat 
a  hub  matrix  as  a  special  kind  of  hub-dominant  matrix. 

Definition  1  A  matrix  A  e  [R',xm  is  called  a  candidate-hub-dominant  matrix  if  its 
first  m—  1  columns  ah  i—  1,. . .,  m—  1,  satisfy  the  following  conditions: 


m—  1 


for  i—  1 1, 


(1) 


where  £>  1  and  17  <  2.  If,  in  addition,  the  last  column  has  its  Euclidean  norm  greater 
than  or  equal  to  that  of  any  other  column,  then  A  is  called  a  hub-dominant  matrix 
and  the  last  column  is  called  the  hub  column. 

For  a  given  hub-dominant  matrix  A,  its  Gram  matrix  is  denoted  by  Q  \—A'A 
and  called  the  system  matrix  associated  with  the  matrix  A.  Assume  that 
A  —  [a\  a2  ■  ■  ■  am_\  am]  is  a  hub-dominant  matrix.  The  corresponding  system  matrix 
Q  is  partitioned  into  a  form  as  follows: 


(2) 
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where  the  matrix  D,  the  vector  c  and  the  real  number  b  in  (2)  are 


'  llaill2 

(aua2) 

(ci\ ,  Clm—  i) 

(d[  5  dm) 

D  = 

{a2,ai) 

1  \a2\\2 

(@2,  dm—\ ) 

,  C  = 

(02, 

-  {dm—  \i  ta\) 

(am—\ ,  af) 

ii  M  2 

1 1  d/fi —  1 1|  _ 

_  (0/77—1  ?  0/77  )  _ 

respectively.  The  matrix  D  represents  the  correlation  among  the  non-hub  nodes. 
The  inequality  (1)  implies  D  being  a  diagonally  dominant  matrix.  This  is  the  reason 
why  A  is  named  a  diagonally  dominant  hub  matrix.  The  vector  c  reflects  the 
correlation  between  the  non-hub  nodes  and  the  hub  node. 

When  the  parameters  £  and  i]  in  (1)  are  identical  and  both  equal  to  1,  we  have 
J2"=i  .(Oi,aj) |  =  ||o,||2,  which  imply  that 

(at,  af)  —  0  for  all  1  <  i  <  j  <  m  —  1 . 

This  means  that  all  non-hub  columns  are  orthogonal  to  each  other  with  respect  to 
the  Euclidean  inner  product.  In  such  case,  the  matrix  D  in  (2)  and  (3)  reduces  to  a 
diagonal  matrix  with  ||ai  ||2,  ||aj  ||2, . . . ,  \\am-\ ||2  as  its  diagonal  entries.  The  matrix  A  is 
called  a  hub  matrix  [1]  and  the  corresponding  Q  is  an  arrowhead  matrix  [3,4]. 

It  is  of  particular  interest  in  estimating  the  eigenvalues  of  the  system  matrix 
associated  with  a  hub-dominant  matrix,  mostly,  the  ratio  of  leading  two  eigenvalues 
of  the  system  matrix  via  the  elements  of  the  hub-dominant  matrix  [1],  To  this  end,  we 
introduce  definitions  of  hub-gaps  for  a  hub-dominant  matrix  and  eigengaps  for  the 
associated  system  matrix. 


Definition  2  (hub-gap  and  eigengap)  Let  A  e  [R''xm  be  a  matrix  with  its  columns 
denoted  by  arranged  in  increasing  order.  For  the  z-th 

hub-gap  of  A  is  defined  to  be 


HG,(d) 


\\a„,-i+i  I 


0/77—7 


Let  Q  be  the  associated  system  matrix  with  eigenvalues  denoted  by  A.1; . . . ,  Xm  with 
*1  <  ■  ■  ■  <  For  i—  1, . . . ,  m  —  1,  the  /-th  eigengap  of  Q  is  defined  to  be 


EG, (0 


bm—i+\ 

l-m—i 


With  this  notation,  an  estimate  of  the  eigengap  of  the  system  matrix  Q  was  given 
in  [1]  in  terms  of  the  hub-gaps  of  the  corresponding  hub-dominant  matrix.  We  state 
the  estimate  in  the  following  theorem. 


Theorem  1  [7]  Let  A  e  J5''xm  be  a  hub-dominant  matrix  with  its  columns  denoted  by 
ai,a2, . .  ■  ,am  arranged  in  increasing  order.  Let  Q  —  A'Ael R"'xm  be  a  corresponding 
system  matrix  with  its  eigenvalues  denoted  by  X i,...,Xm  with  Ai  <  •  •  •  <  Xm.  Let  du  and 
a,  denote  the  diagonal  entry  and  the  sum  of  the  magnitudes  of  off-diagonal  entries, 
respectively,  in  row  i  of  D  for  i  —  1, . . . ,  m  —  1 .  Then 

*  HGiM)  <  EGi(0  <  +  b  +  £'=  i  07 .  (4) 

2  d(m—2)(m—2)  &m—2 

We  make  three  remarks  on  this  theorem. 
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Remark  1  We  would  like  to  take  this  opportunity  to  emphasize  that  the  difference 
d(,n-2)(m-2)  —  &m-2  in  the  denominator  of  the  upper  bound  shown  in  (4)  should 
be  understood  as  the  second  largest  number  in  the  sequence  dn—  au 

dl2  —  ’  “(m—  l)(m— 1)  —  am- 1- 

Remark  2  The  expression  in  (4)  fails  to  give  an  upper  bound  of  EG\(Q)  when 
d(m—2)(m—2)  2  —  t). 

Remark  3  When  A  is  a  hub  matrix,  (4)  becomes 

^HGi(zi)  <  EGj(g)  <  (HG^)  +  1)HG2GQ. 

In  comparison  to  the  lower  and  upper  bounds  of  EG[(Q)  given  by  Theorem  4  in  [1], 
we  see  that  the  upper  bound  of  EG[(Q)  in  [1]  is  the  same  as  the  one  in  the  above 
inequality,  but  the  lower  bound  doubles.  In  this  sense  there  is  room  to  further 
improve  the  lower  bound  of  EGi(Q)  in  (4). 

These  remarks  partially  motivated  us  to  conduct  our  current  work  in  this  article. 
Since  the  eigengap  EGi(g)  is  defined  as  the  ratio  of  leading  two  largest  eigenvalues 
of  Q,  it  is  desirable  to  estimate  the  leading  two  largest  eigenvalues  of  Q  in  terms  of 
the  entries  of  the  associated  hub-dominant  matrix  A  with  a  high  accuracy  in  order 
to  have  tighter  lower  and  upper  bounds  of  EG[(Q).  This  will  be  our  main  task  in  the 
next  section. 


3.  Some  lemmas 

Let  A  be  a  hub-dominant  matrix  and  Q  be  the  associated  Gram  matrix  or  the  system 
matrix.  The  matrix  Q  has  a  form  of  (2),  which  is  partitioned  according  to  the 
non-hub  columns  and  the  hub  column  of  A.  We  will  use  this  structure  of  Q  to 
estimate  its  eigenvalues. 

When  A  is  a  hub  matrix,  Q  is  an  arrowhead  matrix.  The  Cauchy  interlacing 
theorem  is  usually  adopted  to  estimate  the  eigenvalues  of  an  arrowhead  matrix  [5]. 
In  [2],  we  provide  new  estimates  for  the  eigenvalues  of  arrowhead  matrices.  It  was 
proved  there  that  the  estimated  lower  and  upper  bounds  of  the  eigenvalues  are 
tighter  than  that  by  using  the  Cauchy  interlacing  theorem.  In  the  meanwhile,  it  was 
numerically  examined  that  the  estimates  for  the  least  and  largest  eigenvalues  are 
better  than  that  based  on  a  method  in  [6],  Because  of  this,  we  propose  to  apply  the 
results  in  [2]  for  new  estimates  of  the  eigenvalues  of  Q  when  A  is  a  hub-dominant 
matrix. 

The  approach  we  adopt  is  motivated  from  the  structure  of  Q  in  (2).  By  looking 
at  it,  when  the  submatrix  D  in  Q  is  diagonal,  Q  is  an  arrowhead  matrix.  Therefore 
the  results  in  [2]  can  be  directly  applicable  in  such  situation.  Fortunately,  if 
the  submatrix  D  is  not  a  diagonal,  D  is  unitarily  similar  to  a  diagonal  matrix  F,  i.e. 
there  exists  an  orthogonal  matrix  U  such  that 

D  =  U'YU, 

where  Y  —  diag(/x ! ,  /x2, . . . ,  lim-\)  with  /U-i  <  fi2  <  •  *  •  <  Hence,  we  have 


T  z' 

’ U 

~D  c~ 

- 

_z‘  h  _ 

1. 

_  c‘  b  _ 

1. 
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where  2  =  Uc.  Therefore,  the  eigenvalues  of  the  matrix  Q  are  the  same  as  the 
arrowhead  matrix 


B  := 


fr 


(5) 


We  call  B  an  arrowhead  matrix  induced  from  the  system  matrix  Q. 

In  what  follows,  the  eigenvalues  X ,■  of  a  symmetric  matrix  i,elR"x"  are  always 
ordered  such  that  A.i  <  X2  <  ■  ■  •  <  X„.  When  the  dependence  of  the  eigenvalues  on  the 
matrix  P  needs  to  be  determined,  we  simply  write  A](P)  <  X2(P)  <•••<  X,,(P)  instead. 

Next,  we  present  the  bounds  for  the  eigenvalues  of  the  arrowhead  matrix  B 
developed  in  [2]. 


Lemma  1  [1 1]  Let  Be  I 


and 


be  an  arrowhead  matrix  in  (5)  having 


(6) 


m— 1)  With 

Ml  —  P'2  — 

—  Mm— 

1  <  b.  Then 

min  { /r  1 ,/-  (yyf) } 

if  ;=  i; 

bj(B)< 

Pj 

if  2  <  j  <  m  — 

AM 

if  j  —  m 

f-(f) 

if  .7  =  1; 

Xj(B)  > 

maxj/x;_ 

•f-m 

if  2  <  j  <  m  — 

/+(¥) 

if  j  —  m, 

(V) 


where 


f±(p )  =  pb  +  ^l-p)b±  J(l-p)2b2  +  4\\z\A . 

As  a  direct  application  of  Lemma  1 ,  we  have  the  following  result  on  the  bounds 
of  the  eigengap  of  the  system  matrix  Q. 

Theorem  2  Let  A  e  [R"xm  be  a  hub-dominant  matrix  with  its  columns  denoted  by 
a\,...,am  arranged  in  increasing  order.  Let  Q  —  A'AeM"'xm  in  (2)  be  the 
corresponding  system  matrix  and  B  in  (5)  be  the  arrowhead  matrix  induced  from  Q. 
Then  we  have 

f+(p i/b)  .  f+Ohn-i/b) 

Mm— 1  niax{/Xm_2,y-(Mm-l  /^)} 

where  /xi,/x2, . . . ,  /rm_i  arranged  in  increasing  order  are  the  diagonal  elements  of  the 
diagonal  matrix  T  in  B. 

Corollary  1  Let  A  e  IR"xm  be  a  hub-dominant  matrix  with  its  columns  denoted  by 
ai,...,am  and  arranged  in  increasing  order.  Suppose  that  all  non-hub  columns  a\, 
a2,  ■  ■  ■ ,  am_\  are  orthogonal  to  each  other  with  respect  to  the  Euclidean  inner  product. 
Then  the  lower  and  upper  bounds  ofEGfQ)  in  Theorem  2  are  sharper  than  those  in 
Theorem  1. 
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Proof  When  all  non-hub  columns  a\,  a2, . . . ,  am_\  are  orthogonal  to  each  other  with 
respect  to  the  Euclidean  inner  product,  we  know  that 

AO  =  du  —  \\an\\2,  for  i—  1, . . . ,  m  -  1. 


In  this  case, 

{b  +  dn)  +  J{b-dn)2  +  4\\c\\2  b 

- — - SI - -  j 

that  is,  the  lower  bound  of  EGi(0  in  Theorem  2  is  sharper  than  that  in  Theorem  1. 

To  show  the  upper  bound  of  EGi(0  in  Theorem  2  is  sharper  than  that  in 
Theorem  1,  i.e. 

l)(m—  \)/b)  ^  t/(m— l)(m— 1)  T  b 

max) d(m—  i )(m—  I  )■]  -  (d(m—  1 )(m—  1) /b ) }  <f(m—  1  )(m—  1 ) 

we  need  to  show  that 


f  +(d(m—\)(m—\) l b)  <  d(m—  l)(m— 1)  T  b.  (9) 

It  is  equivalent  to  showing 

ll^ll  ^  d(fn-l)(m—l)b- 

The  above  inequality  is  true  because 

m—  1  m—  1  ii  m2 

l|c||2  =  ^2,\{aham)\2  =  {a;-,am>  |2 

i=l  ;=  1  Ill'll 

m—  1  | 

^  ll^m— ill  ^  TT  (#/?  Clm)  I 

tr  n^'ii 

^  1 1  @m—  ill  II  1 1  —  l)b. 

This  completes  the  proof.  ■ 

From  Corollary  1,  we  have  tighter  bounds  for  the  eigengap  EG[(Q)  than  that  in 
Theorem  1  when  A  is  a  hub  matrix.  For  A  being  a  hub-dominant  matrix,  the  bounds 
given  in  Theorem  2  involve  the  smallest  eigenvalue  ii\  and  two  largest  eigenvalues 
p;„_2  and  /rm-i  °f  D  which  are  not  available  in  general.  Hence,  we  need  to  estimate 
these  eigenvalues  in  terms  of  entries  of  the  associated  hub-dominant  matrix  A 
explicitly.  To  this  end,  for  the  eigenvalues  and  /rm_2,  we  need  to  give  lower  bounds 
of  them;  for  the  eigenvalue  /rm-i,  we  need  both  lower  and  upper  bounds. 

The  following  lemmas  are  useful  in  estimating  eigenvalues  ac,  Atm-2  and  aGj-i 
of  D. 

Lemma  2  [4,  p.  191]  Let  PeK"x"  be  a  symmetric  matrix,  let  r  be  an  integer  with 
1  <  r  <  n,  and  let  Pr  denote  any  r  x  r  principal  submatrix  of  P  ( obtaining  by  deleting 
n  —  r  rows  and  the  corresponding  columns  from  P).  For  each  integer  k  such  that 
1  <k  <r  we  have 


h(P)  <  h :(Pr)  <  Xk+n-,{P). 
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This  lemma  shows  that  the  eigenvalues  of  P  can  be  estimated  through  eigenvalues 
of  its  principal  submatrices.  To  make  this  lemma  useful,  the  key  is  to  identify 
principal  submatrices  of  P  whose  eigenvalues  are  easily  computed.  Clearly,  we  can 
explicitly  calculate  the  eigenvalues  of  the  principal  submatrices  Pr  for  r—  1  and  2. 
This  leads  to  the  following  result. 

Lemma  3  Let  P  e  1R',X"  be  a  positive  semi-definite  symmetric  matrix.  Then 

(1)  M(P)  <  mini<;<„ pa  and  X„{P)  >  max,  <,<„/?„, 

(2)  k„_i (P)  >  maxi<(<i<„  |  ((pa  +  pjj)  -  pu  -  pjf)1  +  4 pr.), 

(3)  X„(P)  >  max,<(</<„ \ (( Pii  +  pj})  +  ,J(Pii  ~  Pjj)2  +  4P2)- 

Proof  The  first  item  is  a  well-known  result.  It  is  a  direct  consequence  of  Lemma  2 
by  setting  r—  1.  To  prove  the  results  in  the  second  and  third  items,  we  set  r  —  2  in 
Lemma  2.  The  corresponding  matrix  P2  has  the  following  form: 


L  aH  aii  J 

where  i  and  j  are  two  distinct  integers  between  1  and  n.  Notice  that  all  diagonal 
elements  of  P  are  non-negative.  The  smallest  and  largest  eigenvalues  of  P2  are 

\  Upu  +  Pii)  ~  y/(P»~Pjj)2  +  4Pij)  and  \  (( Pu  +  Pii )  +  y/iPa+Pii)2  +  44') . 

respectively.  The  results  stated  in  items  2  and  3  are  followed  from  Lemma  2.  ■ 

Since  \{{pu  +  pjj)  +  ^{pn  -  pjj)1  +  4 />?•)  >  ma ApiuPjj},  the  lower  bound  for  a„(P) 
from  item  3  of  Lemma  3  is  tighter  than  that  from  item  1 . 

Lemma  4  Let  P  e  IR"X'!  be  a  diagonally  dominant  and  positive  semi-definite 
symmetric  matrix.  Let  us  denote  er,-  the  sum  of  the  magnitudes  of  the  off-diagonal 
elements  of  the  i-th  row  of  P,  i.e.  a,  —  Ylj=ij^i  \Py\-  us  define 

X  :=  diag(pn  -  a\,p22  -a2,...  ,pnn  -  <?„)■ 

Then 


HP)  >  H*) 


(10) 


for  i—  1, 2, . . . ,  n. 

Proof  Since  P  is  a  diagonally  dominant  and  positive  semi-definite  symmetric 
matrix,  then  all  diagonal  entries  of  X  are  non-negative.  Define 

Y:=P-  X. 

Clearly,  Y  is  a  diagonally  dominant  symmetric  matrix  and  its  diagonal  elements  are 
a  j  which  are  non-negative  for  /  =  1,2...../;.  Therefore,  Y  is  a  diagonally  dominant 
and  positive  semi-definite  symmetric  matrix,  hence  all  eigenvalues  of  Y  are 
non-negative.  Since  P  is  a  sum  of  the  diagonal  matrix  X  and  the  non-negative 
symmetric  matrix  Y,  Weyl’s  monotonicity  theorem  (see,  e.g.  [7,  page  63]  or 
[4,  page  182])  implies  (10).  ■ 
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Both  Lemmas  3  and  4  provide  a  lower  bound  for  X„(P).  Since 

ma x[pu  -  ay  i  =  1, 2, . . . ,  n]  <  pnn  <  *  f(pm  +  pjj)  +  yj (pim  -  pjj)2  +  4 p;tj 

for  all  j—  1, . . .  ,n,  we  choose  the  lower  bound  of  A.„(P)  provided  by  Lemma  3. 

Regarding  the  lower  bound  of  based  on  Lemmas  3  and  4,  we  have  the 

following  result. 

Lemma  5  Let  PeM"xn  be  a  symmetric  diagonally  dominant  and  positive  semi- 
definite  matrix.  Let  as  denote  at  the  sum  of  the  magnitudes  of  the  off-diagonal  elements 
of  the  i-th  row  of  P,  i.e.  a ,■  =  Ylj=i  jyi  \Py\-  Then  the  lower  bound  of  A.„_i (P)  given  in 
Lemma  3  is  tighter  than  that  in  Lemma  4. 


Proof  We  assume  that  the  diagonal  entries  pa  of  P  have  been  arranged  in  increasing 
order.  We  consider  two  different  cases.  In  the  first  case,  we  assume  that 

Pnn  -  er„  =  max{  pa  -  or.  z  =  1,2, ... ,  n}, 


then 


max  \  ( (pa  -\-pjj)  -  J {pu  -  Pjj)1  +  4pl  \  >  max  {plt  -  \pa\)  >  max  (p„  -  a,), 

\<i<j<n  Z  \  v  /  l<z<«— 1 

where  max , i ( pu  —  cri)  is  the  second  largest  number  among  pu—cr i, 
P22  —  o2, . . . , Pnn  —  rr„,  therefore,  is  the  exact  lower  bound  of  >.„-\(P)  provided  by 
Lemma  4. 

In  the  second  case,  we  assume  that  for  an  integer  1  <jo<n—  1, 


then 


PM  ~  °>o  =  max{  pu  -ay  i=  1,2,...,  n]. 


max  l-  ( (Pii  +Pjj)  -  J  ( Pu-Pjj )2  +  4pl)  >  max  i 

—  Pjojo  ~  l/’roC/o+hl 

^  Pjojo  -  W)  =  max{/z,7  -  ay  i  =  1 , 2, . . . , «} . 

From  both  cases,  we  conclude  that  the  statement  of  this  lemma  is  true.  ■ 

Lemma  6  Let  PeWxn  be  a  diagonally  dominant  and  positive  semi-definite 
symmetric  matrix.  Then 

n 

K(P)  <  max  y>,|<2  max  \pa\. 

1  <i<nL — i  1  <i<n 

r=i 

Proof  Note  that  all  eigenvalues  of  the  positive  semi-definite  symmetric  matrix  P  are 
real  and  non-negative.  The  first  inequality  is  due  to  Gershgorin’s  disc  theorem  while 
the  second  inequality  is  due  to  P  being  a  diagonally  dominant  matrix.  ■ 

Lemma  7  Let  Q  e  [R'"x"'  be  a  positive  semi-definite  symmetric  matrix.  We  partition 
Q  as  follows: 


C P«  +Pjj)  -  ^(Pa-Pjj)2  +  4P1j 


Q  = 


~D 

c' 


c 

b 
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where  D  e  )  anrf  c  e  [fg"'_1.  Let  us  denote  a,  the  sum  of  the  magnitudes  of  the 

off-diagonal  elements  of  the  i-th  row  of  D.  i.e.  ay  :=  \dij\,  where  dy  is  the  ( i,j)-th 

entry  of  D  (of  Q  as  well).  Assume  that  D  is  a  diagonally  dominant  matrix,  then 

K,(Q)  <b+  max  {, du  +  <x,}. 

1  <i<m—\ 

Proof  As  we  know,  from  the  beginning  of  this  section,  the  matrix  Q  can  be 
converted  to  an  arrowhead  matrix  B  given  in  (5)  by  a  similarity  transformation. 
Hence,  the  eigenvalues  of  Q  and  B  are  the  same.  Further,  by  the  fact  that  the  trace  of 
a  matrix  is  equal  to  the  sum  of  its  eigenvalues,  we  have 

AM  +  AD  +  •  •  •  +  AOh-1  +  b  —  X\(Q)  +  X2  (0  +  •  •  •  +  X„,(Q),  (11) 

where  am,  AG,  •  •  •  ,  AG1-1  arranged  in  increasing  order  are  the  eigenvalues  of  the 
diagonal  matrix  r  in  (5). 

Applying  the  Cauchy  interlacing  theorem  for  B  yields 

HQ)  <  ^  <  Vi(e)  (12) 

for  j—  1,2, ...  ,/m—  1.  Applying  (12)  to  (11)  leads  to 

AM  +  AM  +  ■  ■  ■  +  A Lm-\  +  b  >  X\(Q)  +  AM  +  •  •  •  +  A Lm-2  +  Xm(Q). 


Hence, 


Xm(Q)  <  b  +  Atm-t  —  X\(Q). 

Since  ki(0  >  0  and  D  is  a  diagonally  dominant  matrix,  the  result  of  this  lemma 
follows  from  the  above  inequality  together  with  Lemma  6.  ■ 


4.  Lower  and  upper  bounds  of  the  eigengap  EGx(0 

In  this  section,  we  utilize  the  results  given  in  the  previous  section  to  provide  new 
lower  and  upper  bounds  for  EGi(0,  where  Q  is  the  system  matrix  associated  with  a 
hub-dominant  matrix  A.  We  further  show  that  the  new  bounds  are  tighter  than  that 
provided  by  Theorem  1. 


Theorem  3  Let  AeM!’xm  be  the  hub-dominant  matrix  having  a\,  02,  ■  ■  ■ ,  am  as  its 
columns  with  || <3 1 1|  <  ||a2||  <  ■  •  •  <  ||am|| .  Let  Q  —  A1  A  e  IR'"xm  in  (2)  be  the  correspond¬ 
ing  system  matrix.  Let  us  denote  Oj  l^yl  as  sum  of  the  magnitudes  of  the 

off-diagonal  elements  of  the  i-th  row  of  D  =  [diff'f}\  ■  Then 


U(s/b) 

max  ]  y  ']j  j  |  djj  | 


<  EGK0, 


(13) 


where  s  —  min  {da  —  of.  i—\,...,m—\}. 

Proof  Let  B  in  (5)  be  the  reduced  arrowhead  matrix  of  Q  and  let  am,  AG,  ■  ■  • ,  AG7-1 
be  the  diagonal  elements  of  the  diagonal  matrix  T  in  B  which  have  been  arranged  in 
increasing  order.  From  Theorem  2,  we  have 


EGi(0  > 


f+(d\/b) 

l^m—  1 
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By  Lemma  4,  we  have  /xi  >s.  Since  f+  is  an  increasing  function,  we  get 
f+(ji\!b)  >/+(#)■  By  Lemma  6,  we  have  /im-t  <  max,  J2"=i'  \dy\-  Hence 

f+(i/b) 


EGi  (0  >  > 


Mm-l  max  ,  Eh  141 

This  shows  that  (13)  holds. 

Corollary  2  Under  the  same  conditions  as  in  Theorem  3,  we  have 


Mhg^  +  Jhg^O  ,  J 


4||c|| 


EGi(0). 


(14) 


^{m—  \){m—  1)  / 

Proof  Since  H  e  D?"x”'  is  a  hub-dominant  matrix,  the  submatrix  D  of  Q  in  (2)  is 
diagonally  dominated.  Hence  £>0  and 

f+UJb)  > /+( 0)  =  \  (b  +  >  +  4||,||2)  =  i  f  i  +  >  +  4||c||2) . 

The  fact  ||z||  =  ||c||  is  used  above.  Since  D  is  a  diagonally  dominated  matrix,  we  get 

m—\ 


max 


}  '  \dij\  ^  max  2 da  — 

z f  1</<W—  1 


7=1 


Hence, 


- ./+^//4  >  h+J^  +  4\\c\\-  =  1  /HGi(^}  +  JucffA)  ■ 

maxi<i<m_i  V  =i  4  4 iX,„_ i)  4\  y 


4||c||2  ' 

d(m-l)(m-\)  / 


Because 


1  HC,M+ Ihg:m+4M 


>^HG,(0, 


(ra—  l)(m—  1)  / 

thus,  Theorem  3  provides  a  better  lower  bound  for  EG,(0  than  Theorem  1. 

Next  we  turn  to  give  an  estimate  of  the  upper  bound  of  EG,(0. 

Theorem  4  Let  AeWn*m  he  a  hub-dominant  matrix  having  a\,a2,...,am  as  its 
columns  with  ||(Zi||  <  ||a2||  <  ■  •  •  <  \\am\\ .  Let  Q  =  A1  A  e  IR'"X”'  in  (2)  be  the  correspond¬ 
ing  system  matrix.  Then 

min  {b  +  sx,f+(sx/b)} 


where 


EG'(0  mzx.{s2,f-(sT,/b)} 

m—  1 

■?i  :=  max  V*  \d»\, 

l<i<m—  1  — — ■f 

7=1 

s2  :=  max  ^  ( (da  +  djj)  -  7(4  -  dj)2  +  4 4  ) , 

\<i<j<m—\  Z  \  v  / 

•?3  :  =  max  \( ( 4  +  4)  +  7(4  -  4)2  +  44  )  • 

\<i<j<m—\  Z  \  v  7  / 


(15) 


11 


Downloaded  by  [Lixin  Shen]  at  04:58  08  July  201 1 


Applicable  Analysis 


1237 


Proof  Since  EGi(G)  is  the  ratio  of  km(Q)  and  (Q),  we  need  to  give  an  upper 
bound  for  Xm(Q)  and  an  lower  bound  for  respectively. 

By  Lemma  1,  we  know  km(Q)<f+(nm_i/b).  By  Lemma  6,  we  get  pm-\<S\. 
The  monotonically  increasing  property  of /+  implies 

MG)  </+(*/*)■  (16) 

On  the  other  hand,  by  Lemma  7  we  obtain  another  upper  bound 

K,(Q)<b  +  sx.  (17) 

Hence,  combining  (16)  and  (17)  together  yields 

MG)  <  min{b  +  s\,f+{s\/b)}.  (18) 

Now,  we  develop  the  lower  bound  for  km_\ (Q).  From  Lemmas  3-5,  and  the 
comments  following  Lemma  4,  we  know  that 

Pm- 2  >  s2  and  (im_  i  >  .s'3 .  (19) 

By  Lemma  1  and  the  monotonically  increasing  property  of/_,  we  have 

Mi(G)  >  max{ /r ,„_2 ,/_ (pm~ \/b)}  >  max{s2,f-(s3/b)}.  (20) 

This  completes  the  proof.  ■ 

Remark  4  When  A  is  a  hub  matrix,  Equation  (9)  implies  that  min{ft  +  5l5 
f+(s\/b)}  —f+(s\/b).  Therefore,  the  upper  bound  for  EGi(G)  given  by  Theorem  4 
is  the  same  as  the  one  given  by  Theorem  2. 

Corollary  3  Under  the  same  conditions  as  in  Theorem  4,  the  upper  bound  for 
EGj(G)  provided  by  Theorem  4  is  tighter  than  the  one  given  by  Theorem  1. 

Proof  By  Lemma  5,  we  have 


d(m—2)(m—2)  C7m— 2  L  ■t,2,  (2 1) 

where  07  :=  l<4'l  is  the  sum  °f  the  magnitudes  of  the  off-diagonal  elements  of 

the  7-th  row  of  D  in  (2).  The  difference  u,(m-2)(m-2) —  erm_2  in  above  inequality  should 
be  understood  as  the  second  largest  number  among  du  —  o\,  d22  —  o2,  ■  ■  ■ , 

d(m—  I  )(m—  I )  CJm—  1 . 

Next,  we  will  show  that 


m—2 

^1  ^  d(jn— 1)(/77 —  1)  ^  ^  &j- 

7=1 


(22) 


We  consider  two  different  cases.  In  the  first  case,  we  assume  that  M”  , 1  \d(m-\)j\  —  s\. 
Then 


m—2  m—2  m-2 

—  d(m— l)(m—  1)  +  ^  ']  \d(m—l)j\  —  d(m—l)(m—l)  T  ^  '  dj(m—  \ )  <  d(m—  l)(m_l)  T  ^  CTj. 

7=1  7=1  7=1 
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In  the  second  case,  we  assume  that  there  is  an  integer  i0  between  1  and  m  —  2  such 
that  Jfj=\  \dioj\  =  si .  Hence, 

m—2  m—2  m—2 

■*i  —  di0i0-\-  }  '  |t/jQ/|  +  \dj0(m—\)\  <d(m~ i)(m_i)+  y  '  | djj0 1  +  ay0  <  ^  ' °j ~b • 

7=  1  7=1  7=1 

./-Tb  74  7)  747i 

Therefore,  (22)  holds. 

By  (21)  and  (22),  we  have 

min{h  +  s\,f+(s\/b)}  +  s\  b  +  d(m- +  Yl'jLi  aj 

max{s2,/_(s3/h)}  s2  d(m-2)(m-2)  &m— 2 

This  implies  that  the  upper  bound  for  EGi(g)  given  by  Theorem  4  is  tighter  than 
the  one  given  by  Theorem  2.  ■ 


5.  Hub-dominant  matrices  from  equiangular  tight  frames 

In  this  section,  we  focus  on  designing  a  class  of  hub-dominant  matrices.  We  begin  the 
discussion  from  hub  matrices. 

Let  A  —  [a\  a2  ■  ■  ■  a,„]  be  a  hub  matrix  with  am  as  its  hub  column.  By  the  definition 

of  hub  matrices,  Y17=\  l(a7>a/)l  =  \\a]  II2  f°r  .7=  1 _ ,m—  1,  i.e.  1  in  (1). 

We  ask  whether  we  can  construct  hub-dominant  matrices  satisfying  the  following 
properties: 


m —  1 

^2  \  {Oj,at)\  =  fi\\aj  ||2,  for  j  =  1 . m-  1 ,  (23) 

;=1 

where  1  <  /J  <  2.  The  answer  to  this  question  is  yes.  We  can  construct  a  class  of  such 
hub-dominant  matrices  from  equiangular  tight  frames.  The  definition  of  equiangular 
tight  frames  is  given  as  follows. 

Definition  3  The  set  of  unit  vectors  {v1;  v2, . . . ,  v^},  where  v,  elR',  is  called  an 
equiangular  tight  frame  (ETF)  if  it  satisfies  two  conditions:  (i)  For  some 
nonnegative  a,  we  have  \(vk,  vf  |  =  a  when  1  <k<j <N  and  (ii)  Y^k=ivkvk  —  ^j Id, 
when  Id  is  the  identity  matrix  in  [R'vx'\ 

We  remark  that  the  first  condition  indicates  the  vectors  in  the  set  being 
equiangular  while  the  second  condition  indicates  the  vectors  generating  a  tight 
frame,  that  is  (see,  e.g.  [8]), 


v  =  vk)vk  for  all  ve  OS'.  (24) 

^  k=  1 

Equiangular  tight  frames  have  been  studied  recently  and  have  found  applications  in 
communications  and  coding  theory  (see,  e.g.  [9-12]). 

Immediately,  the  definitions  of  hub-dominant  matrix  and  ETF  lead  to  the 
following  results. 
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Proposition  1  Fix  m  >  3  and  let  A  —  [a\  a2  •  ■  ■  an,]  be  a  hub-dominant  matrix,  where 
a,e  IR"  for  1  <  /  <  m.  If  the  collection  of  non-hub  columns  ci\,  a2, . . . ,  am_\  of  A  forms 
an  ETF,  then  (i)  ||a,„||  >  1;  (ii)  fi  —  (m  —  2)a  +  1,  where  ft  is  the  constant  in  (23)  and 
a  —  |  (a,-,  af  \  for  all  1  <  i<j  <m  —  1  and  (iii)  a  <  — L. 

We  are  particularly  interested  in  the  case  where  yS  =  2,  i.e.  a  —  m]_2,  under  the 
conditions  of  Proposition  1.  Mercedes-Benz  systems,  a  special  type  of  ETFs,  satisfy 
this  additional  requirement. 

An  ETF  {b\,  b2,...,  bn+  J  in  IR"  is  called  a  Mercedes-Benz  system  if 

{bk,  bj)  —  for  kfj. 

Lemma  8  Let  {b\,  bj,-.-,  bn+ 1}  in  IR"  be  a  Mercedes-Benz  system.  Then 


77+ 1 

/=1 

Proof  For  every  1  <  j  <n  +  1 


n+l 


n+\ 


n+l 


bi ’  bi )  =  (bj>  bj)  +  ( XI  bi,  bj )  =  1  +  (bn  bj)  =  1  -  -  =  0. 


i  i=  1 


[  i=  1 

¥j 


7=1 

¥j 


By  the  second  condition  in  the  definition  of  an  ETF,  we  have 


77+1 


77+1  /  77+ 1 


T=0- 


This  completes  the  proof.  ■ 

We  remark  that  the  concrete  constructions  of  Mercedes-Benz  systems  are 
discussed  in  [13]. 

Lemma  9  Fix  m  >  3  and  let  A  —  [a\  a2  -  ■  ■  am\  e  be  a  hub-dominant  matrix 

with  am  as  its  hub  column.  Let  Q  in  (2)  be  the  system  matrix  associated  with  A.  If  the 
set  of  the  non-hub  columns  a\,a2,...,  am_\  of  A  forms  a  Mercedes-Benz  system,  then 
Z)e[R(,n-|»x(„,-i);  a  suj)mairix  of  Q  in  (2),  has  a  form  of 


„  m— 1,  1  , 

D  — - -Id - -ee , 

m  —  2  /77  —  2 


(25) 


where  e  :=[1, 1, . . . ,  l](e  IR'”  .  Further,  the  eigenvalues  of  D  are  0  f simple )  and 
with  multiplicity  m  —  2. 


Proof  Since  the  set  of  the  non-hub  columns  a\,a2,  ■  ■  ■  ,am-i  of  A  forms  a 
Mercedes-Benz  system,  then  +•,  the  (/,  y)-th  entry  of  D  is 


dij  =  {at,  af 


i 

l 

777  —2 


if  i=j; 
if  ifj- 


It  leads  to  (25). 
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Since  ee ’  is  a  rank-1  matrix  and  has  'jjff  as  its  only  non-zero  eigenvalue  with  e 

as  the  corresponding  eigenvector,  hence  the  eigenvalues  of  D  are  0  (simple)  and 
with  multiplicity  m  —  2.  ■ 

Theorem  5  Fix  m  >  3  and  let  A  —  [ai  a2  -  •  •  am]  e  [fg(,”-2)xm  be  a  hub-dominant  matrix 
with  am  as  its  hub  column.  Let  Q  in  (2)  be  the  system  matrix  associated  with  A.  If  the 
set  of  the  non-hub  columns  a\ ,  go,  . . . ,  am_\  of  A  forms  a  Mercedes-Benz  system,  then 
the  eigenvalues  of  Q  are  0  with  multiplicity  2,  h’/T/z  multiplicity  m  —  3,  and  b  +  jff 
(simple). 

Proof  By  Lemma  9,  there  exists  an  orthogonal  matrix  U e  [R>('”_Ox(m-i>  tjiat 

D  —  U‘T U,  where  T  =  diagto,  — — * , — 'V 

\  m  —  2  m  —  2) 


Many  orthogonal  matrices  U  can  serve  this  purpose.  Here  we  simply  choose  U  to  be 
the  (m  —  1)  x  (m  —  1)  DCT-II  matrix  whose  first  row  is  |  e  (see  [14]).  With  the  help 
of  U,  we  have 


"r  z  ~ 

'U 

D  c~ 

'If 

_z‘  b_ 

1 

c‘  b 

1 

where  b—\\am\\2,  c‘—[(am, a{)  (am,  a2)  ■  ■  ■  (am,  am_  1)]  and  z—Uc.  Write  z  — 
[z  1  z2-  ■  -zm_  1].  Then,  by  Lemma  8, 


zi 


m—  1 


--e  c  = 


m  ■ 


Vi 


m 


-j-  'y  ]  fm  -  < 
1  !=1 


v 


m 


I  V  \ 

(  am,  2_,  ai) 


0. 


By  (24), 


o')  VYl  1 

ll^ll  =  II^CII  =  llcll  = - «l 

m  —  2 

Collecting  all  these  results  together,  we  get 


11  J  = 

m  —  2 


(26) 


det( Aid 


T  z 


=  A  A 


=  A  I  A 


m  —  1 
m  —  2 
m  —  1 
m  —  2 


n—2 


(A  -  b)  —  A  A 


m  —  1 
m  —  2 


n— 3 


m- 3 


A  —  b 


m  —  1 
m  —  2 


This  indicates  that  0  with  multiplicity  2,  with  multiplicity  m  —  3  and  b  + 
(simple)  are  the  eigenvalues  of  Q.  This  completes  the  proof.  ■ 

Under  the  conditions  of  Theorem  5,  we  have 


EG!(0  =  1  +  — — Th¬ 
in  —  1 


For  a  fixed  b  (i.e.  a  fixed  hub  column),  EGfiQ)  increases  when  the  number  of 
non-hub  columns  increase  and  EGi(0  is  bounded  by  1  +b.  For  a  fixed  number 
of  non-hub  columns,  EGi(0  is  a  linear  function  of  b  with  the  slope  ffj- 

From  Theorem  1,  the  lower  bound  of  EGi(0  is  \b.  However,  the  upper  bound 
cannot  be  applied  because  d(, „-2)(m-2)  —  rrm_2  =  0. 
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By  a  simple  calculation,  we  know  that 


5  =  0,  5i 


2, 


S2 


m  —  3 
m  —  2  ’ 


As 


m  —  1 
m  —  2 5 


where  s,  si,  s2  and  ,v3,  are  defined  in  Theorems  3  and  4.  By  using  (26)  and  the 
definition  of  /±,  we  have 


A(sJ  =  i 


4  (in  —  1) 
m  —  2 


f+(si/b)  >b  +  5i,  f-(s3/b)  =  0. 


Then,  the  lower  and  upper  bounds  for  EGi(Q)  given  by  Theorems  3  and  4  are 
the  following: 


1 

4 


4(/71  —  1 )  ^ 
m  —  2 


<  EGi(0)  < 


m  —  2 
m  —  3 


(2  +  b). 


6.  Conclusion 

We  have  developed  improved  bounds  for  the  eigengaps  of  the  system  matrices 
associated  with  hub-dominant  matrices.  It  would  be  interested  in  studying  the 
eigenvectors  of  the  system  matrices  in  the  future  research.  Extension  of  the  current 
results  to  multi-hub  matrices  (i.e.  two  or  more  hub  columns)  will  be  studied  as  well. 
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